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Abstract 

We study QCD effects in single graviton production at the CERN Large Hadron Collider (LHC) 
in the Randall-Sundrum (RS) Model. We present in detail the complete next-to-leading order 
(NLO) QCD corrections to the inclusive total cross sections. The NLO QCD corrections enhance 
significantly the total cross sections and decrease efficiently the dependence of the total cross 
sections on the factorization and renormalization scales. We also examine the uncertainty of the 
total cross sections due to the parton distribution function (PDF) uncertainties. For the differential 
cross sections on the transverse momentum (qt) of the graviton, within the CSS resummation 
formalism, we resum the logarithmically-enhanced terms at small qt to all orders up to NLO 
logatithmic accuracy. Combined with the fixed order calculations, we give consistent predictions 
for both small qt and large qt- 
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I. INTRODUCTION 



Now search for extra dimensions has been one of the major objects at the LHC, since 
its physical effects can appear at the TeV energy scale. The idea of extra dimensions was 
revived in the 1990 'sjl[ Q, 0, S, 0], which can bring new solutions to the gauge hierarchy 
problem and be used to resolve some problems of the SM such as the origin of the fermion 
masses and their hierarchy. 

So far, there have been various extra dimension models, which can be divided into two 
major classes according to the geometry of the background space-time manifold. The first 
one includes the ADD model and its variants, which extend the dimension of the totoal 
space-time to D = 4 + 6, propose a factorizable metric, and get large size of the extra 
dimensions 1/Mp). In the ADD model, the SM particles live in the usual 4— dimensional 
space-time, while gravity can propagate in the additional (5-dimensional space, which is 
assumed for simplicity to be compactified on the 5-dimensional torus with a common 
radius R. Then the 4-dimensional Planck scale Mp is related to the fundamental scale Mg 
as follows |l[ Q]: 

= Mf+2(27ri?)^ (1) 

where Ms ~ TeV. According to Eq. (Q), deviations from the usual Newtoniangravitational 
force law can be expected at distances smaller than i? ~ 2 x 10^^'' x 10^cmj7|. For 5 > 2, 
ADD is consistent with the current experiments jSj] since gravitational forces are not yet well 
probed at distances less than about a millimeter. However for 6 = 2, there are constraints 



arising from, e.g., supernova cooling, which require Mg > 10^— lOOTeV if 5 = 2 [2|- 

The second one includes the 5-dimensional RS model 2] and its variants, in which a 
warped metric is introduced and the size of the extra dimension needs not to be too large 
compared with the Planck length. In the RS model, the extra dimension is assumed to be 
an S1/Z2 orbifold, which has two fixed points, ^ = and 9 = n. At each fixed point, there 
is a 3-brane, and the brane a.t 6 = n corresponds to the brane we live on, while the one at 
^ = is the high energy brane. Between the two 3-branes is a slice of AdS space, where 
only the graviton can propagate into. Moreover, the 4-dimensional metric is the function of 
the coordinate of the 5th dimension, i.e. 

ds^ = e-''="=l'^l(r/^, + -^h^,)dx>'dx'' - rld(P\ < |0| < tt, (2) 



where is a scale of order of the Planck scale and relates the 5-dimensional Planck scale 
to the cosmological constant, is the compactification radius, and h^y is the graviton. 

AftCT solving the 5-dimensional Einstein equation, we can get the tensions of the two 
branes 7] 

K) = = 12A;M3, (3) 

and from Eq.Q, we can get the relation between and 4-dimensional reduced Planck 
scale MpQ 

M^ = i^(l-e-2^'^='^), (4) 

from which we can see that for moderately large values of the compactification scale rc, the 
relation between M* and Mp almost does not depend on rc, and it is completely different 
from the results in the ADD model. Compared with the ADD model, the RS model present 
a different solution to the gauge Hierarchy problem: the physical mass m of a field on the 
brane where our world live on, is related to the fundamental mass parameter itiq as following 

m = e"''''"''mo, (5) 

thus the hierarchy problem can be solved if kvc ~ 12. 

In the RS model, there also exist KK towers of massive spin-2 gravitons which can interact 
with the SM fields, and we have the following 4-dimensional effective Lagrangianjo, 10 1: 

1 1 oo 

^ = -=^T"^(a:)/.S(x) - -T'^^ix) ^ C^(^) (6) 



where 

A. = e--"Mp = (7) 

Xik 

and is at the electroweak scale. Thus the coupling of the massless graviton h^^^ is suppressed 
by the Planck scale, and the ones of the massive graviton h^^'> by A^r, but which is only TeV. 
The masses of the nth graviton KK excitation modes are also at the electroweak scale, which 
are given by 

m„ = fcx^e-'^'^^" = mi — , (8) 

Xi 



where the the nth. roots of the first order Bessel function. 



From Eqs.Q and (jH)), the graviton sector of the RS mode^ 



the two parameters mi and k/ Mp. Current constraints 



IC 
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is completely determined by 



12| | for the parameters of the 



RS model are from the theoretical requirement, the low energy precise measurement and 
also the data from Tevatron, from which we expect 0.01 < k/Mp < 0.1 and < lOTeV. 

There are two classes of effects that can be used to probe extra dimension in the RS model 
at high energy colliders: real graviton emission and virtual KK tower exchange. In the RS 
model, the lightest massive graviton can have a mass of several hundred GeV, and may 
be produced copiously at the LHC. More importantly, it has much larger couplings to the 
SM particles than the ones in the ADD model, thus it may decay into observable particles 
and hence be detected. And there have been detailed analysis [13| which demonstrate that 
using channels pp — > h^"^^ e~^e~, 77..., we can probe the massive graviton in the RS model 
with masses up to several TeV. However, those analysis [l^l are based on the LO results, in 
order to improve the precision of the theoretical predictions, the higher order QCD effects 
are necessary. In the ADD model and the RS model, the NLO QCD corrections to the 
virtual graviton production at the LHC have been discussed in Ref. 2, however, the K 
factors contributed from different parts, the scale dependence and the PDF uncertainty for 
above processes needs further studies. Moreover, they also did not consider the kinematic 
distribution of the events, which is very important in designing the strategy of discovery. In 
this paper, we study the transverse momentum distribution of the massive graviton at NLO 
in QCD, and all order soft gluon resummation effects on the distribution to give reasonable 
predictions. 

The arrangement of this paper is as follows. In Sect. II, we show the LO results and 
define the notations. In Sect. HI, we present the details of the calculations of both the 
virtual and real parts of the NLO QCD corrections. In Sect. IV, we give the transverse 
momentum distribution. In Sect. V, we present the detailed numerical results for the total 
cross sections and also the transverse momentum distribution. Sec. VI contains a brief 
conclusion. 



II. LEADING ORDER CALCULATIONS 



The related Feynman diagrams which contribute to the LO amplitude of the partonic 
process gp{pi)gtiP2), QriPi)QsiP2) h^uJ are shown in Fig. ^ The relevant Feynman rules 



can be read easily from the ones in the ADD model presented in Ref. 



1^ . from which 
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FIG. 1: Leading order Feynman diagrams for pp ^ h 



in) 



we can get the LO amplitude in n = 4 — 2e dimensions as following 
<' = X 



Pi ■ P2C, 



pi/, pa I pu,p(T I ^pu,pa 
,A—n 



e-Me'MeZiPi + P2) 



V{P2 



4A. 

lp{Piu - P2u) + lu{pip - P2p) - 2?7^^(^i- ^2 - 2mg) 



(9) 



^pi)<:(pi+P2), (10) 



where 5ab and 5rs are color tensors (a, 6 are the color indices of the initial state gluons, and 
r, s are the color indices of the initial state quarks), C^^^p^, D^^^p^ and Ep_^p^{pi,p2) are the 
coefficients in the couplings between the graviton and gluons, which can be found in Ref . Q| , 
Hr is a mass parameter introduced to keep the couplings dimension 
For the polarization sum of the massive graviton, we havejl^. Il4| 



ess. 



s=l 



where 



(11) 



PpuafB = ^ [vpc^VuP + VpPVua - ^ ^ ^ VpuVo^P^ + . . . , 

the dots represent terms proportional to the graviton momentum A;^, and since k'^T^ 
give a vanishing contribution to the amplitude. For convenience, below we define 

2 2 

X 



(12) 

= 0, 

(13) 



n-1 3 - 2e 

Moreover, in order to avoid introducing external ghost lines while summing over the gluon 
helicities, we limit ourselves to the sum over the physical polarizations of the gluons jisl ]. 



I.e. 



Pr = E4iki)eUh) = -g'-' + 



(14) 



rp Tl'i ' ki ij^i ' kj^ 

where the index i (=1,2) labels the two external gluons, and rii ^ ki is an arbitrary vector. 
This polarization sum obeys the transversality relations 



ki^P^ — kiu — fli/iP^ — P^ f^iu — 0. 



Thus we can get the relevant partonic cross sections as following: 



^Jf^ = ^27r5(.-m^)|M^;; 



(2 



e vr 



32A2 



5(1 -f), 



(15) 



(16) 



(LO) 
99 



-2'n5{s - ml)\M, 



1 -e)7r 



5(1 -f) 



(17) 



25-" V- -n.,-«« , 24A2 

where s = {pi+ ^"^^ ^ = °^nl ^■ 

The LO total cross sections at the LHC are obtained by convoluting the partonic cross 
sections with the parton distribution functions (PDFs) Gq^q^g/p in the proton: 



(LO) = (LO) , {LO) 
w — " gg ' '^gq ' 



a 



(J. 



(LO) 
99 

(LO) 



1 .1 1 
dxi / dx2- 

TQ Jto/xi 2 

1 .1 

dxi dx2 

TO Jtq/xi 



Gg/p{Xl, fJ'f)Gg/p{X2, fif) + Gg/p{X2, fif)Gg/p{Xl, fif 

Gq/p{Xl, fIf)Gq/p{X2, /i/) + Gq/p{X2, fif)Gq/p{Xl, 

where Tq = m^/S'o, 5*0 = (14TeV)^ and is the factorization scale. 



a 



(LO) 

99 ' 



cr. 



(0) 
gq 5 



(18) 
(19) 



III. NEXT-TO-LEADING ORDER CALCULATIONS 

The NLO QCD corrections consist of the following contributions: the virtual corrections 
arising from loop diagrams of colored particles, the real contributions arising from the radi- 
ation of a real gluon or a massless (anti) quark, and the contributions of mass factorization. 
In the following, we will calculate these contributions separately. We use dimensional regu- 
larization (DREG) in d = 4 — 2e dimensions to regulate the ultraviolet (UV) and infrared 
(IR) divergences. 

For the partonic cross section, the total virtual corrections can be written as 

±.y ■^L.unren i " con ^90^ 

99,qq ~ 99m ~'~ 99,qq^ \ ' 

6 



where the first part in the right hand contains the radiative corrections from the one-loop 
vertex diagrams, and the second part is the contributions from the counterterms involving 
only the the wavefunction renormalization constant for the external fields. 

The 0{as) virtual corrections to the partonic cross section can be expressed as 

- - ^) 

X _ , ^ 1 2 _ 203 357v\ 

U e?^ 48 eiR + 24 e,^ + 8 96 + 288 > ^ > 



a 



247rA2 ^ ' \ 6 ejji 3 eiR 18 6 ^ ^ 



with Ce = rxi^2e) (^^^^ ' which are UV finite, but still contain the IR divergences. Here, 
the IR divergences include the soft divergences and the coUinear divergences. The soft diver- 
gences are canceled after adding the real emission corrections, and the remaining collinear 
divergences can be absorbed into the redefinition of PDF 17j , which will be discussed in the 
following subsections. 

The real corrections consist of the contributions from the radiation of a real gluon or 
a massless (anti) quark. The Feynman diagrams for the real gluon emission sub-process 
gg,qq — ^ h^^J are shown in Figs. |21 and Fig. IHl respectively. The Feynman diagrams for 
massless (anti) quark emission (the diagrams for the antiquark emission sub-processes are 
similar and omitted here) are shown in Fig. 0] 

For the real gluon emission sub-processes g{pi)g{p2), <i{Pi)q{P2) —* g{P3)h^f2)^ ^^e partonic 
cross sections is 



'•real 
99m 



I J iM-'/rfF^, (23) 



with 



+e\tu'^[26 - 9(4 - 2e) - 2(-5 + 2(4 - 2e))] - ~^K ^[16 - 6(4 - 2e) + 4e\] 
~ ^^-t^[16 - 6(4 - 2e) + 4e^x] - '^^ ~^^ g^(t + u)[16 - 6(4 - 2e) + 4e\] 



' ^ ~ ^\^[16 - 6(4 - 2e) + Ae^x] - ~ (^^ + tu + u'^^lW - 6(4 - 2e) + 4e\] 



2 L ^ ' 4s 



+|[2etM(26 - 9(4 - 2e) - 6ex) " ^(1 - e)(t' + u^)[l6 - 6(4 - 2e) + 4e\]}, (24) 




FIG. 2: Feynman diagrams of real gluon emission sub-processes gg hfj,J + g. 




FIG. 3: Feynman diagrams of real gluon emission sub-processes qq — > hnJ + g. 

-(e - 3)[(e - l)t'^ + 2(3e - 2)ut + (e - l)u^]s'^ 

-{t + u) [(e - 1) V + (6e(e - 3) + 8)ut + (e - 1) V]s 

-2tu[{e - 1)^2 + 26-^^ + (e - [e(ex + 3) - 2] I , 
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(n) 

FIG. 4: Feynman diagrams of real quark emission sub-processes gq — > 



where \M^^°-^^\ is the squared matrix of the real gluon emission sub-processes, in which the 
colors and spins of the outgoing particles have been summed, and the colors and spins of 
the incoming ones have been averaged over, and the final state 2-body phase space is 

dT, = — -i (^)\mi - ff X ^-(1 - v^dv, (26) 

87rr(l - e) V / 



where 



y = i(l + cos^), (27) 
t={p,-p^f = -s{l-f){l-v), (28) 
u={p2-p^f = -s{l-f)v, (29) 

and 9 is the the angle between the incoming gluon and the outgoing gluon. 

Combining the contributions of the virtual corrections and the real gluon emission, we 
still have the collinear divergences, which can be absorbed into the redefinition of the PDF 
at NLO, in general called mass factorization [17]. This procedure in practice means that 
first we convolute the partonic cross section with the bare PDF Ga/p{x), and then rewrite 
Ga/p{x) in terms of the renormalized PDF Ga/p{x, ^jyvn the numerical calculations. In the 
MS scheme, the scale dependent PDF Ga/p{x, jif) is given bvjist 

Ga/pix,fif) = Ga/pix) ^f71fr^(~~^) ] / ^Pai3iz)Gi3/p{x/z), (30) 



/3 



where Pa/^iz) are the leading order AltareUi-Parisi sphtting functions jioj ]. 

Here, for the real gluon emission sub-processes gg h^^J and qq h^f^J, we first consider 
only the contributions from Pgg and Pgg, and we can get the relevant counterterm arising 
from the PDF redefinition as following: 



5a,, = 2x^C:%^.p(J)(.), (32) 



where _ p 



27r ' 24A2 

^(1-2!) (^) and z = m^/(xiX25') = f. 
Summing up the virtual, real emission and PDF redefinition contributions, we have the 
IR finite results 

=^55 +^gg + '^^gg= 32A2 



203 35^ .h>(l-f)N 

^ 12 36 ^ ^ ' V 1-f A 



1-f , /(l-f)2\ Inf 3 11 3 llf2 , , 

+6[_l + _ + t(l_f)„„(L__L)_6_____ + _tH._f. (33) 

-TVLO ±.real , ±.V , x;^ (l^^)o^s^"^n 

c^OT = + + -^^^g- = 24A2 ^^T" 



X 



4 /m2\ l + f2 3^,. 4 TT^ 16/ln(l-f) 



4,^ (1-tY 8 Inf 16 16f2 

--1 + r In- 



3 ' ' f 3 1 - f 9f 9 

For the real (anti-)quark emission sub-processes gq{q) Q{Q)h^^J, the relevant results 
can can be got in the similar way as above. First, the partonic cross sections is 



'•real 



I J \M-^/d^2, (34) 



with 



= -^|M^-f 'I {s ^ t), (35) 



72 



l^sT'l = -~T^jM',t\ ^ (36) 
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where jiW^^'^y is the squared matrix of the real (anti-)quark emission sub-processes, in 
which the colors and spins of the outgoing particles have been summed, and the colors and 
spins of the incoming ones have been averaged over. 

Secondly, the relevant counterterm arising from the PDF redefinition for real (anti-)quark 
emission sub-processes are shown as following, where we consider only the contributions from 

Pgg and Pgg-. 

Summing up the virtual, real emission and PDF redefinition contributions, again we get 
the IR finite results for real (anti-)quark emission sub-processes. 



;-NLO( ;-NLO\ 



9Q 



96A2 



a— 



x|[4 
+ [4 



l + (l-f) 



T 



mi 



2 , ,2m i„M / \ _j 

2 f 2 



+ ((l-f)2 + f2)]ln 



(38) 



r V r 

The NLO total cross section for pp /i^") in the MS factorization scheme is obtained by 
summing up the Born, virtual, real emission and PDF redefinition contributions. In terms 
of the above notations, we have 



a 



[NLO) 



(LO) ^ JLO) ^ ^{NLO) ^ ^{NLO) ^ ^(NLO) ^ JNLO) 



(LO) , (LO) 



qg 



gq 



gq 



+ ( dXl I dX2-Gg/p{Xl,IJ,f)Gg/p{X2,fif)+Gg/p{X2,IJ,f)Gg/p{Xl,fif) 

+ f dXl j dX2 Gg/p{Xl, Hf)Gg/p{X2, fXf) + Gg/p{X2, iJ,f)Gg/p{xi, Hf) 

Jto Jto/xi I 

+ f dXl I dX2 Gg/p{Xl, flf)Gg/p{X2, flf) + Gg/p{X2, flf)Gg/p{Xl, Hf) 

Jto Jtq/xi I 

+ f dxi f dx2 Gg/p{xi,fif)Gg/p{x2,Hf) + Gg/p{x2,Hf)Gq/p{xi,Hf) 

Jto Jto/xi I 



a 



(NLO) 



gg 



{NLO) 

qq 



.(NLO) 

gq 

.(NLO) 

gq 



(39) 



Finally, we note that our NLO results in Eqs.(|H^ and are the same as the ones in 
Eq.(3.38) of the second paper in 

Ref.Q, 

except the differences by overall factors. 
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IV. TRANSVERSE MOMENTUM DISTRIBUTION 



In this section we investigate the transverse momentum distribution of the massive gravi- 
ton. At LO the graviton is kept at zero qt due to momentum conservation and the dis- 
tribution is proportional to 5'^{qT)- Thus the LO distribution at non-zero Qt belongs to 
where momentum conservation is retained by the additional parton emitted. The 
distribution can be obtained from the squared amplitudes of the real emission processes, i.e. 
Eqs. (j211), and (jH^j) . However, the corresponding fixed order result of the transverse 
momentum distribution is only valid when qx is not too small compared with the mass of the 
massive graviton If qj- <^ m„, the corresponding parton emitted would be either soft 
or collinear to one of the initial partons. Thus, large logarithms like ln(m^/g|,) will appear 
and will dominate over the cross section for sufficiently small qx- In general, there should 
be double logarithms for each gluon attached to the initial quarks due to the overlap of soft 
region and collinear region. As a result, the perturbative expansion would be controlled by 
asln^(m^/gy) rather than as- The convergence of the perturbation series will be spoiled if 
asln^(m^/gy) approaches unity. In order to make use of the perturbation theory with the 
existence of large logarithms at each order, one must reorganize the perturbative expansion 
to resum the large terms. In this paper, we use the Collins-Soper-Sterman (CSS) resum- 
mation formalism [2^ to calculate all order soft gluon effects on the transverse momentum 
distribution. 

In the CSS formalism, the differential cross section we are considering can be written as 



da 



(total) 



da 



(resum) + F(gT, a;?, a;^). 



(40) 



dq^dy 



dq^dy 



where 



da 



(resum) = ^ 

a,f3 




(resum) 



(41) 



dqj'dy 




(42) 



ab 



and the resummed part can be expressed as an inverse Fourier transformation 




(resum) 




/ 



d'^bexp (ib ■ qx] Wapih., m, x?, X2) 




(43) 
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with 



X exp 



/i2 



(44) 



where a/3 = ^f^f, qq, ah = gg, qq, q{q)g, b is the impact parameter conjugating to qr, Jo is 
zero order Bessel function of the first kind, and = e^m„/A/i, a;^ = e'^rrin/ y/s. Here 
Ci{i = 1, 2, 3) are constants of order 1 which are by convention j20] chosen to be 



Ci = Cs = 2e'^^ = bo, C2 = 1, 
and / is the convolution of the PDFs and the coefficient functions C 

fa/hix,fi) =J2 I —Ca^iz,asilJ,))f^/hix,fi), 

and the coefficients A, B and C can be expanded to series in 

^1 Vvry 



(45) 



(46) 



n=0 



(47) 
(48) 
(49) 



and they can be calculated order by order in perturbative theory. In our case, since the 
massive graviton is colorless, thus the lowest order coefficients is the same as the ones in the 
case of gg ^] and Drell- Yan 22| • For the gg channel, we have 



A(^) = 2N^ = 6, 5^^) = -2po = (33 - 2nj)/6, 



and for the qq channel, we have 



5« = -'-c, = -2, 



Ci:'>{z) = 6^p5{l-z). 



(50) 



(51) 
(52) 



With these coefficients, we can actually sum up all terms like a^L^""^ and a^L^""^. 

However, the resummed part is still not able to be calculated perturbatively. The reason 
is that in Eq. the integral over the impact parameter b extends to infinity, while the 
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integrand involves the strong coupling constant and the PDFs at scale h^/h^ where they 
are not well defined if h is large enough so that ho/h enters non-perturbative region. Collins, 
Soper and Sterman, in their original paper [2^, suggested that one can use a cut-off 6max 
and regard the effects from h > femax as non-perturbative input. Practically, they replaced 
W{b) in Eq. (gSD by 

Wib) = Wih)FMb), (53) 

where 

h 

(54) 



max J 



and -Fnp(6) parameterizes the non-perturbative effects. Since 6* never exceeds 6max, W^(&*) 
can be calculated perturbatively, and the theoretical uncertainty mainly relies on the func- 
tion -Fnp. Recently, Landry, Brock, Nadolsky and Yuan (BLNY) |23] proposed the form 



NP 



exp < — 6^ 



^ ■ on, 



fi'i + fi'2 In 7^ + 9\9z ln(100a;ia;2 



(55) 



They take 6max = 0.5GeV ^, Qo = l-6GeV and the parameters giii = 1,2,3) are fitted to 
the available Drell-Yan data, which are given by 

51 = 0.21, 52 = 0.68, 53 = -0.60. (56) 

The another term in Eq. (PT|) . the Y term, is the remaining contributions which are not 
resummed. Since it contains no large logarithms, it can be reliably computed by fixed order 
truncation of the perturbative series 

where the first term in the right hand is the fixed-order perturbative results, and the second 
term is the asymptotic part of the differential cross section, defined as the terms which are 
at least as singular as l/q^ when ^ 0, which can be got by expanding the resummed 
part, i.e. Eq. pHj) . In our case, we have 

^(asym) = ^<^|{/,/p(x?, /.,)/,/p(x°, ;.,) (2iV.ln^ - 2/?o) 

+ (^99 ° f)g/pixi^ f^f)fg/pix2, /^/) + fg/piA f^f)iPgg ° f)g/pixl, f^f) 

+ (x?^x°)j, (58) 
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+ {Pqq O f)q/p{X2,fif)fq/p{xl,fif) + f q/ p{xl, li f){Pqq O f)g/p{x1,fif) 



+ {x1^xl)\, (59) 



^^"^^ (asym) = fg/piA ^^f){Pqg ° f)g/pi4^ f^f) 



^/p(x?,/i/)(Pgg o /)g/p(a;°,^/) 



1 a^7o_ 



+ (60) 



where 0"°^ = ^^5- and cr°^ = g^r- The result of ^^^{cisym) is similar to the one of 
-^pj^{asym), and thus omitted here. 



V. NUMERICAL RESULTS 

As mentioned in Section I, there are two additional free inputs in the RS model: mi and 
k/Mp, which have the following constraints: 0.01 < k/Mp < 0.1 and = ^liMp < 10 TeV. 
In our numerical calculations, for convenience, we choose the input parameters as A^r and 
nil. For A^ = 4 (8)TeV^, from current constraints, we have 150 GeV < rrii < 1.5 TeV 
(300 GeV < mi < 3 TeV). 

Moreover, for the NLO total cross sections a'-^^*^^ and the contributions from different 
parts (including (tJ^^^\ 4?"^°^ and (xj^^^) + a^,^^^^), the NLO (MS) PDFs [2^ is used 
throughout this paper. For the LO results, we define two cross sections as following: 

^(Loi) . j^Q paj-^onic cross section convoluted with NLO (MS) PDFs; (61) 
^(L02) . j^Q pg^j-^QYiiQ cross section convoluted with LO PDFs, (62) 

and correspondingly two K factors: 

O-(NLO) (j(NLO) 
^1 = ^(LOiy' ^2 = ^^Lo^. (63) 

As the above definitions, Ki measures only the size of the NLO QCD corrections to the 
cross sections, while K2 accounts for the effects of changing parton distribution functions 
additionally. As for the renormalization and factorization scales, we always choose /i^ = 
fif = rrin, unless otherwise specified. 
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In Fig. El we show the dependence of the total cross sections of pp h'^^} at the LHC 
as functions of mi, assuming = 4TeV. The NLO and LO total cross sections decrease 
when mi increases. The LO total cross sections are in general over several pb, and reach 
100 pb when mi = 500 GeV. Moreover, the figure also shows that the NLO QCD corrections 
enhance significantly the LO total cross sections, which are in general several tens percent. 

Fig. ini shows the dependence of the K factors on mi, based on the results in Fig. El When 
mi varies from 500 GeV to 1.5 TeV, the Ki factor ranges from 1.46 to 1.44, and the K2 
factor ranges from 1.61 to 1.71. In addition, we give the different parts of the Ki factor, 
which show that the contributions from cr^^^^^ ranges from 0.51 to 0.44, the contributions 
from (Tqf^'^'' range from 0.03 to 0.14, and the contributions from (^g^^^^ ~^'^'g^^^^ range from 
-0.08 to -0.14, respectively. 

Fig. 13 shows the dependence of the total cross sections of pp — > h^^} at the LHC on the 
renormalization scale (/i^) and the factorization scale (/i/), assuming /i^ = /i/, A^r = 4 TeV 
and mi = 1 TeV. The scale dependence of the NLO total cross section is much smaller than 
that of the LO cross section. For example, the LO cross sections 

~ ±17.8% (~ ±17.3%), when /x^ = /ij ranges between 500 GeV and 4 TeV, while the NLO 
ones vary by ~ ±9.3%. Moreover, we also give the scale dependence of the different parts 
in a^^^'^\ for example, when Hr{= l^f) ranges between 500 GeV and 4 TeV, (^^^^'^^ varies 
from 6.9pb to 10.4pb, CTq^^'^^ varies from 0.56pb to 0.93pb, and cr^^^*^^ + '^5?^^*^'' varies from 
-0.5pb to -3.5pb, respectively. 

To estimate the uncertainties in the total cross sections due to the uncertainty of PDFs, 
we take the 41 sets of CTEQ6.1 PDFs to calculate the LO and NLO rates hA. Fig.Hshows 

n 

the PDF uncertainties (defined here as the Eq. (3) in Ref. [2^) in the LO and NLO total cross 
sections for pp h^^} production at the LHC, as functions of mi, assuming = 8 TeV. It 
turns out that the PDF uncertainties in the LO and NLO total cross sections increases as 
mi increases. Moreover, when mi is small (< 1.5 TeV), the PDF uncertainties in the LO and 
NLO total cross sections are about the same, while when mi becomes large (> 1.5 TeV), the 
NLO rate has a larger uncertainty than the LO rate due to the PDF uncertainties, especially 
at large mi. 

Figs, inland El gives the transverse momentum distribution of h^^} frompp — >■ h^i^} process 
at the LHC, assuming A^ = 4 TeV, for mi = 1 TeV and 2 TeV, respectively. The peaks of 
the distribution appear at about 18 GeV and 13 GeV. The differential cross sections decrease 
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sharply with the increase of gr, which indicates that most events will happen in the relatively 
low qt region, where the resummation effects are essential. Moreover, we also plot the various 
parts of the differential cross sections in Eq.(|4Up. The perturbative and the asymptotic cross 
sections agree very well at small transverse momentum. On the other hand, the resummed 
and the asymptotic part are not canceled completely at high due to the higher order 
effects included in the resummed one, so that the total one and the perturbative one will 
differ at large qt- This can be considered as the theoretical uncertainties. In principle, one 
can return to the perturbative result for qx > g™*, where g™* is arbitrarily chosen in the 
intermediate gr region. However, in order to make the transition smooth, one must introduce 
some kinds of matching procedure which could also lead to uncertainties. In our work, as 
shown by Eq. (j57|) . we subtract from the Y term the expansion of the resummed part of the 
same perturbative order, and this matching procedure between small qT and large qT region 
prevents double-counting of perturbative results and also guarantees a uniform theoretical 
accuracy over the entire qr region j^]. 

From Eqs. ()42|1 and (jKTj) . we know that the transverse momentum distribution can 

be divided into three parts, i.e. the contributions from gg, qq and gq + gq channels: 

"^"^ -(total) = (total) + (total) + ^^^f±^ (total) (64) 



dqj^dy dqj.dy dqj^dy dqj-dy 



where 



dq^dy ^^'^^^^^ dq^dy^^^^^^ ~'~ dq^dy^^^^^^ 

In Figs imiT^ andlT^ we thus plot these three contributions to the the transverse momentum 
distribution of h^^^ based on the results in FigsEl respectively. We can see that for all these 
three parts, the perturbative and the asymptotic cross sections agree very well at small 
transverse momentum. 
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VI. CONCLUSIONS 



In summary, we have calculated the next-to-leading order total cross section and trans- 
verse momentum distribution of single massive graviton production at the LHC in the RS 
model, including all-order soft gluon resummation effects. Our results show that the LO 
total cross sections are in general over several pb in most of the parameter space, and can 
reach 100 pb when mi = 500 GeV. The NLO corrections enhance significantly the total cross 
sections, which is in general several tens percent, and reduce efficiently the dependence of 
the total cross sections on the renormalization/factorization scale. We have also examined 
the uncertainty in total cross sections due to the PDF uncertainties, and found that the 
uncertainty in NLO cross sections is shghtly larger than that in LO ones, especially at large 
mi. For the transverse momentum distribution, within the CSS resummation formalism, we 
resum the logarithmically-enhanced terms at small qt to all orders up to NLO logarithmic 
accuracy. Combined with the fixed order calculations, we give consistent predictions for 
both small qt and large qr- Our results can be useful to the simulation of the events in the 
future collider experiments. 

Acknowledgment s 

We thank Yang Gao for useful discussion. This work was supported in part by the 
National Natural Science Foundation of China, under grants No. 10421503 and No. 10575001, 
the Key Grant Project of Chinese Ministry of Education under grant No. 305001. 



[1] N. Arkani-Hamed et al., Phys. Lett. B429, 263 (1998); N. Arkani-Hamed et ak, Phys. Rev. 

D59, 086004 (1999); N. Arkani-Hamed et ak, Phys. Lett. B436, 257 (1998). 
[2] L. Randall, R.Sundrum, Phys. Rev. Lett. 83, 3370 (1999); L. Randall, R.Sundrum, Phys. 

Rev. Lett. 83, 4690 (1999). 
[3] J. Lykken, Phys. ReV. D54, R3693 (1996). 
[4] E. Witten, Nucl. Phys. B471, 135 (1996). 

[5] P. Horava and E. Witten, Nucl. Phys. B460, 506 (1996); Nucl. Phys. B475 (1996) 94. 
[6] I. Antoniadis, Phys. Lett. B246, 377 (1990). 



18 



[7] C. Csaki, TASI Lectures on Extra Dimensions and Branes, "hep-ph / 0404096'. 

[8] E. G. Adelberger [Eot-Wash Group Collaboration] liep-ex/0202008 

[9] J.L. Hewett, Phys. Rev. Lett. 82 4765 (1999). 

[10] H. Davoudiasl, et al., Phys.Rev.Lett. 84, 2080 (2000). 

[11] JoAnne Hewett, Maria Spiropulu, Ann.Rev.Nucl.Part.Sci. 52, 397 (2002). 

[12] Prakash Mathews et al., JHEP 0510, 031 (2005); Prakash Mathews et al., Nucl.Phys. 
B713, 333 (2005); Prakash Mathews, V. Ravindran, |hep-ph/0507250l M. C. Kumar et al., 



hep-ph/0604135" 



[13] B. C. Allanach, et al., JHEP 0212, 039 (2002). 

[14] T. Han et al., Phys. Rev. D59, 105006 (1999). 

[15] W. Beenakker et al., Phys. Rev. D40, 54 (1989). 

[16] G. 't Hooft, M. J. G. Veltman, Nucl. Phys. B44, 189 (1972). 

[17] G. Altarein et al., Nucl. Phys. B157, 461 (1979); J. C. Colhns et al., in: Perturbative Quantum 

Chromodynamics, ed. A. H. Mueller (World Scientific, 1989). 

[18] B. W. Harris, J. F. Owens, Phys. Rev. D65, 094032 (2002). 

[19] G. AltarelU, G. Parisi, Nucl. Phys. B126, 298 (1977). 

[20] J. C. Collins, D. E. Soper and G. Sterman, Nucl. Phys. B250, 199 (1985). 

[21] R.P. Kauffman, Phys. Rev. D44, 1415 (1991). 

[22] P.B. Arnold, R.P. Kauffman, Nucl. Phys. B349, 381 (1991); C. Balazs et al., Phys. Lett. 

B355, 548 (1995). 

[23] F. Landry et al., Phys. Rev. D67, 073016 (2003). 

[24] J. Pumplin et al., JHEP 0207, 012 (2002). 

[25] D. Stump et al., JHEP 0310, 046 (2003). 

[26] J. Pumplin et al., JHEP 0207, 012(2002). 

[27] G. Bozzi, B. Fuks and M. Klasen, Phys. Rev.D74, 015001 (2006). 



19 




FIG. 5: Dependence of the total cross sections for the first KK graviton excitation mode direct 
production at the LHC on mi. 



20 





1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 

(NLO), (L02) 

a la 


(NLO)/ (L01) 

^ A =4TpV 

l\ — 1 C V 

71 - 

(NLO), (L01) 

a la 

99 


r (NLO) (NLO),, (L01) 

L"^„„ "^^^TT (NLO), (L01) 

gq gq a _ la 

' ■ 








1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 



500 600 700 800 900 1000 1100 1200 1300 1400 1500 

m. (GeV) 



FIG. 6: Dependence of the K-iacioi on mi, based on the results in FigO 



21 




FIG. 7: Dependence of the total cross sections for the first KK graviton excitation mode direct 
production at the LHC on fir = fJ-f- 
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FIG. 8: The PDF dependence of the total cross sections for the first KK graviton excitation mode 
direct production at the LHC, as functions of mi. 
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FIG. 9: The transverse momentum distribution of the first KK graviton excitation mode from 
PP fijlu process at the LHC, assuming A^r = 4TeV and mi = 1 TeV. 



24 



0.01 



1E-3 



i 1 i 1 1 1 i 1 1 1 i 1 1 


1 ' 1 ' 1 ' 




total . ~\ 




. pert - 


resum 






asym 



1 — I — 1 — I — 1 — I — 1 — I — 1 — I — 1 — I — 1 — I — 1 — I — 1 — I — 1 — 

20 40 60 80 100 120 140 160 180 200 

q^(GeV) 



FIG. 10: The transverse momentum distribution of the first KK graviton excitation mode from 
pp — h^fiu process at the LHC, assuming A-j^ = 4TeV and mi = 2TeV. 
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FIG. 11: The gg part of the transverse momentum distribution of the first KK graviton excitation 
mode, assuming = 4 TeV and mi = 1 TeV. 



25 




FIG. 12: The qq part of the transverse momentum distribution of the first KK graviton excitation 
mode, assuming A^r = 4 TeV and mi = 1 TeV. 




FIG. 13: The gq and gq part of the transverse momentum distribution of the first KK graviton 
excitation mode, assuming = 4 TeV and mi = 1 TeV. 
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